Enhancing molecular conversion efficiency by a magnetic field pulse sequence 
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We propose a strategy to enhance tiie atom-to-molecule conversion efficiency near a Feshbach 
resonance. Based on the mean-field approximation, we derive the fixed point solutions of the classical 
Hamiltonian. Rabi oscillation between the atomic and molecular states around fixed point solutions 
and its oscillation period are discussed. By designing a sequence of magnetic field pulses in analogy 
with Ramsey experiments, we show that a much higher atom-to-molecule conversion efficiency can 
be accessed by tuning the pulse durations appropriately. 
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I. INTRODUCTION 

The study of iiltracold molecular gases has absorbed 
much attention since it is bringing the emergence of novel 
physical phenomena and is motivating potential appli- 
cations of special advantages. For instance, the former 
includes phenomenon of BCS-BEC crossover [ij and ul- 
tracold chemistry [3 , while the latter involves precision 
spectroscopy [EfllSl, and the test of fundamental sym- 
metries 0,13 etc. 

To investigate various novel physical features of ul- 
tracold molecules, one needs to convert cold atoms into 
molecules. In recent experiments, the technique of Fes- 
hbach resonance originally introduced to manipulate 
the interaction between atoms plays an important role 
in the creation of molecules. The route of tuning the 
magnetic field can not only affect the atom-to-molecule 
conversion efficiency, but also make the system exhibit 
different fascinating phenomena. To date for most ex- 
periments, one sweeps the magnetic field sl owly across 
a Feshbach resonance to convert fermionic [l^ [HI or 
bosonic [H, [H, [H, [3 atoms into molecules, but 
this always induces heating and additional particle loss 
to the system. To avoid such a problem, one can add a 
small sinusoidal oscillation to a constant magnetic field 
which is far away from the Feshbach resonance [13, [3 ■ 
As long as the modulation frequency closely matches the 
molecular binding energy, the atoms can be converted 
into molecules with less particle loss. For the conver- 
sion of degenerate fermionic atoms into molecules, one 
can also hold the magnetic field near the Feshbach res- 
onance in the regime with positive scattering length for 
several seconds such that a weakly bound molecular state 
emerges [l^. This technique relies on the long lifetime 
of molecules and works only for fermionic atoms. 

Besides the above magnetoassociation approaches, 
atoms can also be converted into molecules through pho- 
toassociation [2^, [2lj, [13, [23l or the stimulated Raman 
adiabatic passage technique [24 [25j . Motivated by a re- 
cent experiment [26| which was originally aimed to in- 
vestigate the coherence property of the atom-molecule 
mixture, in the present paper we propose another strat- 
egy to efficiently enhance the atom-to-molecule conver- 



sion efficiency through applying a sequence of magnetic 
field pulses. The paper is organized as follows. In Sec.lTTl 
we introduce a simplified two-level model describing the 
atom-molecule conversion system. In the mean-field ap- 
proximation, we derive the classical Hamiltonian and the 
corresponding evolution equations. In Sec. IIIIl the fixed 
point solutions of the classical Hamiltonian are presented. 
Rabi oscillation between the atomic and molecular states 
is discussed, and the dependence of the oscillation pe- 
riod on the system parameters is given analytically. In 
Sec. IIVI we simulate the evolution of the system under a 
sequence of magnetic field pulses in analogy with Ram- 
sey experiments. Based on the discussion, we propose a 
strategy to enhance the atom-to-molecule conversion effi- 
ciency by tuning the pulse durations. The summary and 
discussion are given in Sec. |Vl 



II. MODELING IN PHASE SPACE 

To study the atom-molecule conversion near a Fesh- 
bach resonance, we consider the following Hamiltonian, 



+eaV'iV'a + em^L^m + ll'i/'m^a'i/'a + V'm'i/'iV'i), (1) 

where operators ipa ("01) and -0111 (V'm) annihilate (create) 
an atom and a molecule, respectively; parameters C/a and 
refer to the atomic and molecular interactions while 
[/am refers to atom-molecule interaction. Here and em 
denote the energies of atomic and molecular states, and g 
measures the Feshbach coupling strength between atoms 
and molecules. 

Since the total number of atoms and their compounds 
(molecules) is sufficiently large, we can apply the mean- 
field approximation, i.e., {tpa) — \fn\J pa(i)e*^''^*^ and 
(''/'m) = v^V^PmC^ye'^'"^*-' with n being the mean atomic 
density. Here the atomic and molecular populations, pa 
and pni, as well as the corresponding phases, 6'a and 
0m 7 are introduced. The particle- number conservation, 
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|(V^a)P + 2|(-0m)P = n, requires pa + 2pm = 1- Corre- 
spondingly, in the mean-field approximation, the Heisen- 
berg equations of motion for annihilation operators V'a 
and V'm are equivalent to the following evolution equa- 
tions for the phase and population. 



(2ea - Em + 2L/a - C/am) + (4C/am ~ ~ Un,)pm 

(2) 



9 n^Pm ~ 1 

'- COS 6^ . 

2 y/p^ 



Pm = -5(1 - 2pni)VPmSin( 



(3) 



where we introduced a relative phase — 6^ ~ 20a 
and several new notations, — U^n/h, C/m = Uaji/Ti, 

C/am = Ua.mn/fi, g = g^/u/h, ea = Sa/h, and em = Cm/^ 

for simplicity. Note that these newly defined quantities 
all have the dimension of frequency, i.e., Hz. We need 
to mention that Eqs. ^ and ^ were derived under the 
preassumption Pa 7^ and p^ ^ 0, which leaves out 
one possible fixed point solution that will be retrieved in 

secnni 

With the help of the canonical conjugate relations. 
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the classical Hamiltonian describing the energy of system 
is obtained 



-(2ea - Cn 
4(4C/am 



+ 2C/a - U, 
4t/a - U^)p 



am J Pm 
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+g{l - 2p,ji)^/p^cos9. 



(5) 



For completeness, we included the interaction terms in 
our model Hamiltonian ([T]). However, in most experi- 
ments, these interactions can be typically ignored as long 
as the magnetic field is not so close to the Feshbach res- 
onance. In the following, we neglect those interaction 
terms for simplicity and choose e^. = without loss of 
generality. Additionally, we take < so as to favor 
the production of molecules. 



III. FIXED POINT SOLUTIONS AND RABI 
OSCILLATIONS 

We can now derive the fixed point solutions of the clas- 
sical Hamiltonian ([5]) which is related to Rabi oscillations. 
To determine the fixed point solutions, we set 6 — and 
Pm = as we did in Ref. [131, and obtain 



9 a^Pn 
■ cost' 



1 



2 y/Pm 

- 2pm)^/fhrism0 



0, 
0. 



(6) 



Solving Eqs. (O gives 6* = or tt for 7^ or 1/2. Each 
case corresponds to two possible solutions. The other so- 
lution left out from Eqs. ^ is p^ — 1/2 (the value of 9 is 



not well defined) which is obtained by directly solving the 
original evolution equations for (ipg) and {ipm}- As afore- 
mentioned, Eqs. ([ni do not yield such a solution because 
the assumption p^ 7^ or 1/2 (i.e., pa = 0) in deriv- 
ing Eqs. (O excluded this case. Considering the physical 
constraint < pm < 1 /2 , we summarize in Table U all 
the physical fixed point solutions. 

TABLE I: Fixed Point Solutions. 
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Em > -V^g 


Undefined 


1/2 


Always 



We are interested in Rabi oscillations initiated from a 
pure atomic state, i.e., Pm(0) = 0. Note that this initial 
condition is accompanied with the fact that the value 
of 9 is not well defined at t = 0. It is clear that the 
time-evolution of the system, as shown in Fig. [ija), is 
characterized by the oscillation between the atomic and 
molecular states. 
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FIG. 1: Time-dependence of Rabi oscillations for (a) em/g = 
—0.8 and (b) em/g = 0. (c) The period of Rabi oscillation ver- 
sus 6m /g, where em < is assumed. Note that all quantities 
have already been renormalized to be dimensionless. 

The period of Rabi oscillation can be easily evaluated 
by making use of the conservation of the classical energy 
([5]) whose value is determined to be zero by the initial 
condition. We obtain the following analytical expression 

f 99 

jRabi = <P iTTTrMPm (7) 
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which corresponds to the maximum of the reachable 
molecular population a-lso the nonzero divergent 
point of dO/dHc- The dependence of Tftabi on the value 
of era/ 9 is plotted in Fig. [TJc), from which we can find 
that the period of Rabi oscillation trends to diverge as 
Icm/sl diminishes to zero. As shown in Fig.[T]Jb), no oscil- 
lation is observed at the resonance em = 0. Instead, the 
evolution brings the system from the pure atomic state 
to the molecular one completely. 



IV. MOLECULAR CONVERSION UNDER A 
SEQUENCE OF MAGNETIC FIELD PULSES 

Motivated by Ramsey experiments which are oriented 
to probe the coherence properties between atoms and 
molecules, here we propose a strategy to enhance atom- 
to- molecule conversion efficiency by a specially designed 
magnetic pulse sequence. As shown in Ref. [2g|, the pa- 
rameter e„i can be tuned by varying the time-dependent 
magnetic field B{t) due to em(i) = ^J'dis[B{t) — Bo] with 
Bq being the position of the exact resonance and /idiff the 
difference between the magnetic moments of the atomic 
and molecular states. We suggest tuning the value of 
in the manner shown in Fig. [2fb). This corresponds to 
the idealized pulse sequence consisting of three periods 
of constant magnetic field. Ti and T3 are the durations 
of the first and the second pulses, respectively, and T2 
is the evolution time between the two pulses. Note that 
we have already replaced the ramps of the first and the 
second pulses by abrupt changes of the magnetic field for 
theoretical convenience. 

Let us discuss the evolution of the system initiated 
from the pure atomic state. As shown in Fig. O the 
two magnetic field pulses correspond to e„i{t)/g = —0.8 
when < i < Ti and (Ti + T2) <t<{Tx+T2+ T3). 
From Fig. [T]^a) and [TJc) we can find that, the Rabi os- 
cillation period is glRabi — 3.7 if only the first magnetic 
field pulse exists. Thus in Fig. [2] we choose Ti = TRabi/2 
in order to reach the maximum of the molecular conver- 
sion at the end of the first magnetic field pulse, which 
is about ~ 0.27. During Ti < t < {Ti + T2) the 
smallness of em(i) is accompanied with a fast oscillation 
of small amplitude during the evolution. In our simu- 
lation (7T2 = 1.21 during which the system has already 
oscillated for dozens of periods. This specific value of 
T2 is chosen to make sure that, when t > (Ti + T2) the 
system undergoes a slow Rabi oscillation of large am- 
plitude (obeying energy conservation) whose maximum 
value can get close to the limit of a complete conversion, 
i.e., pra = 1/2. This would never happen if one just had 
applied one magnetic field pulse, which is illustrated in 
Fig. [TJa). In our simulation the maximum conversion 
efficiency locates at gtc — 6.25 with Pm{tc) — 0.49. 

Now we interpret this evolution behavior with the help 
of the energy contour of Eq. ^ in the phase space of pm 
and 6. Figure. [3ja) corresponds to e^/g = —0.8 while 
[3]Jb) corresponds to e^n/g = —50.0. The different values 
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FIG. 2: (a) Time-dependence of the atomic and molecular 
populations, (b) Schematic time-dependence of em{t) under 
the sequence of magnetic field pulses. The horizontal line in 
(a) indicates the limit of the complete conversion of molecules, 
i.e., pm ~ 1/2. The two vertical lines separate the evolution 
process into three regions (I, II, and III) corresponding to 
the time-dependent magnetic field pulses. Durations of each 
region is denoted by Ti (i = 1, 2, and 3). 



of Cm yield different energy contour patterns. Here the 
energy conservation plays the major role. 

During the first magnetic field pulse (I), < t < Ti, 
the system oscillates in the Rabi type obeying the en- 
ergy conservation law. In our simulation, the undefined 
6(0) at the beginning of the evolution is chosen to be 
zero as shown in Figure. [3{a). At t = Ti = rRabi/2, 
the sudden change of magnetic field switches the energy 
landscape described by the contour pattern in Fig. [3I^a) 
to that in Fig. ^h) while keeping the instantaneous lo- 
cation in the phase space unchanged (a' corresponding 
to a) at the moment. Note that the system energy is 
changed during this switch. In contrast, during the se- 
quential magnetic field (II) which is far away from the 
resonance, energy conservation still holds, forcing the sys- 
tem to evolve with a fast oscillation of small amplitude. 
Suppose at t = (ri-l-T2) it evolves to the location (3' given 
in Fig. EKb) . When the magnetic field is switched back to 
the original value of the first pulse at the same moment, 
the system locates at a new position /3 (corresponding to 
(3') which is no longer inside the original trajectory of the 
Rabi oscillation induced by the first pulse. Consequently, 
the system has entered another trajectory with constant 
energy. Furthermore, the energy conservation will def- 
initely bring the system close to a complete conversion 
into molecules (III). Note that the evolution direction in 
Figure. [3] can be determined by analyzing Eqs. ^ and 
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FIG. 3: (Color online) Energy contour for (a) em/g = —0.8 
and (b) em/g = —50 in the phase space of pm and 6. The 
arrows show the evolution direction. The extreme points 
(marked by the plus sign) correspond to the fixed point solu- 
tions discussed in Sec. IIIII pm = 1/2 is unmarked due to the 
not well-defined 9. Note that a (/3) and a' (/?') refer to the 
same point in the phase space. 



V. SUMMARY AND DISCUSSION 

We considered the atom-to-molecule conversion near a 
Feshbach resonance and proposed a strategy to enhance 
the conversion efficiency. In the mean-field approxima- 
tion, we derived the classical Hamiltonian and the cor- 
responding evolution equations, together with the fixed 
point solutions of the classical system. Starting from 



the pure atomic condensate, we investigated the Rabi os- 
cillation and presented the analytical expression for the 
oscillation period versus em/ 9- In analogy with Ram- 
sey experiments, we applied to the system a sequence of 
magnetic field pulses. We found that, by tuning the du- 
rations of each pulse and the evolution time elaborately, 
it is possible to approach a complete conversion that can 
never be achieved by a single pulse. With the help of 
the energy contour in the phase space of and 9, we 
gave a simple interpretation of such evolution based on 
the energy conservation law. 

The advantage of our proposal lies in the fact that it 
does not require the magnetic field to cross or closely 
approach the Feshbach resonance, which can avoid the 
substantial system heating and particle loss induced by 
the enhanced interactions. In our proposal the magnetic 
field in each pulse (we denote it by Bp) could take val- 
ues away from Bq. Thus, we can achieve much higher 
ultimate phase space densities (nA^j^ with Ath being the 
thermal de Broglie wavelength) and hence higher possible 
conversion efficiencies. In our simulation, the magnitudes 
Bp of the two pulses are taken to be the same, but they 
can also be different to obtain a higher conversion effi- 
ciency. However, as \Bp — Bo\ increases, the value of 
decreases accordingly, which may lead to the suppression 
of the amplitude of Rabi oscillation during the two mag- 
netic field pulses. The value of Bp should not to be too 
far away from Bq. For example, the analysis of Fig. O 
together with the fixed point solutions in Table |T] suggest 
that the optimal range of -Bp can be chosen such that 
-V2 < ejg < 0. 

The fast ramps of magnetic field pulses may result in 
the loss of atoms and the production of an additional 
component which is considered as an atomic burst con- 
sisting of correlated pairs with a comparatively high rel- 
ative velocity [1^. To take this effect into account, our 
model should be extended to include the atomic burst 
as a new component of continuum spectrum. This may 
bring about the damping in the Rabi oscillation j28j] , but 
the effects of the atomic burst would decrease as the tem- 
perature drops close to zero. 

The work is supported by NSFC Grant No. 10674117, 
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